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Abstract: Gauge/gravity duality has proved to be a very successful tool for describing strongly 
coupled systems in particle physics and heavy ion physics. The application of the gauge/gravity 
duality to quantum matter is a promising candidate to explain questions concerning non-zero tem- 
perature dynamics and transport coefficients. To a large extent, the success of applications of 
gauge/gravity duality to the quark-gluon plasma is founded on the derivation of a universal result, 
the famous ratio of shear viscosity and entropy density. As a base for applications to condensed 
matter physics, it is highly desirable to have a similar universal relation in this context as well. A 
candidate for such a universal law is given by Homes ' law: High T c superconductors, as well as some 
conventional superconductors, exhibit a universal scaling relation between the superfluid density at 
zero temperature and the conductivity at the critical temperature times the critical temperature 
itself. In this work we describe progress in employing the models of holographic superconductors to 
realize Homes' law and to find a universal relation governing strongly correlated quantum matter. 
We calculate diffusive processes, including the backreaction of the gravitational matter fields on the 
geometry. We consider both holographic s-wave and p-wave superconductors. We show that a par- 
ticular form of Homes' law holds in the absence of backreaction. Moreover, we suggest further steps 
to be taken for holographically realizing Homes' law more generally in the presence of backreaction. 

Keywords: Holography and Condensed Matter Physics (AdS/CMT), 
Gauge-Gravity Correspondence. 
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1. Introduction 

Gauge/gravity duality has proved to be a valuable tool for exploring strongly coupled regimes of 
field theories. The best studied example so far for applications to experimentally accessible strongly 
coupled systems is the application to the quark-gluon plasma. A very important example for this is 
the derivation of the famous result for the ratio of the shear viscosity and the entropy density [1], 

S 4-7T fc B 

Here the physical constants H and fee are written out explicitly in order to illustrate the influence 
of quantum mechanics and thermal physics. 1 It has been shown [2-4] that this result applies 



1 Subsequently, we will set the physical constants c, h and fee to one in the following sections except for 
Section 2.3 and Section 2.1. 
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universally for any isotropic gauge/gravity duality model based on an Einstein-Hilbert action on 
the gravity side. Exceptions are found by considering higher curvature corrections [5] or anisotropic 
configurations, see [6,7] and [8]. Recently, the focus of applying the tools of the gauge/gravity duality 
has been widened to other strongly coupled systems in physics, especially to problems in condensed 
matter physics. In particular, significant progress has been made in describing holographic fermions 
(see [9-11] and references therein), superconductors/superfluids (for instance [12-16] and references 
therein) and to some extent also to lattices [17, 18]. For obtaining a solid general framework 
for condensed matter applications of the gauge/gravity duality, it would be very useful to derive a 
universal relation, similar in importance to (1.1), designed in particular for applications in condensed 
matter physics. Interestingly, the result (1.1) may be understood in the context of condensed matter 
physics by a time scale argument. Here, the properties of quantum critical regions [19,20] give rise 
to a universal lower bound 

k B T 

sometimes called "Planckian dissipation" [21] which can be compared to the possible lower bound 
for v/ s given in (1.1). This seems to imply that the "strange metal phase" is a nearly perfect fluid 
without a quasi-particle description as is the quark-gluon plasma, since both cases do not allow for 
long-lived excitations compared to the energy 

- « e, 1.3 

T 

but rather describe a regime characterized by 

n 



T 



C 



(1.4) 



In the case of the quark-gluon plasma, a possible characteristic time scale can be defined by 



r\ ~ tr. (1-5) 

In typical condensed matter problems at quantum critical points, the relevant energy scale e is set 
by the thermal energy e ~ k^T. An interesting, yet unresolved problem, are the high temperature 
superconductors and their possible relation to quantum critical regions. A very interesting uni- 
versality shown by almost all types of superconductors is Homes ' law. As explained in Section 2 
in detail, this shows some connections to quantum critical regions and the "Planckian dissipation" 
time Tfr. Thus, it is an exciting candidate to find a universal relation for strongly coupled condensed 
matter systems [22] where the usual quasi-particle picture seems to fail. 



In this paper, we outline how Homes' law may be implemented in holography. We follow a simple 
approach which allows to demonstrate the validity of Homes' law in the absence of backreaction 
of the matter fields on the geometry. We also find that our straightforward approach requires 
modifications in the presence of backreaction. We present and discuss possible generalizations and 
indicate directions for further research. The paper is organized as follows. In Section 2 we give a 
self-contained exposition of Homes' law and discuss related condensed matter concepts and their 
holographic realization. In the Sections 3 and 4 we consider holographic s-wave and p-wave super- 
conductors, respectively. In particular, we focus on the effects that arise from the backreaction of 
the gravitational interaction with the matter fields onto the geometry, governing the gravity dual. 
For this purpose, we numerically determine the phase diagram where we use the strength of the 
backreaction as parameter in addition to the ratio of temperature and chemical potential. This sets 
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the ground for the calculations of various diffusion constants, which are discussed for the s-wave 
superconductor in Section 3.5 and for the p-wave superconductor in Section 4.3. We focus on the 
critical diffusion and associated time scales, i.e. the diffusion at the critical ratio of temperature 
and chemical potential depending on the backreaction at which the system transits to the super- 
conducting phase. We find that a particular version of Homes' law is satisfied if the backreaction 
is absent, while further work is required for the case with backreaction, as we explain. Finally, in 
Section 5 we summarize our results and give some possible explanations why Homes' law is not 
confirmed in the approach considered here once the backreaction on the geometry is turned on. 
Some extensions of our original setup to remedy this are discussed in Section 6 along with some 
new perspectives on holographic superfluids/superconductors that might be of interest to pursue 
on their own. 

2. Homes' Law 

With the discovery of high temperature superconductors, the new era heralded by the discovery of 
novel phases of (quantum) matter boosted the need for a new understanding of the classification of 
condensed matter systems. The experimental progress in controlling strongly correlated electronic 
systems and the exploration of strongly coupled fermionic/bosonic systems with the help of ultra- 
cold gases presented a new picture of nature that shook the old foundations of traditional condensed 
matter theory, i.e. Landau's Fermi liquid theory and transitions between different phases classified 
by their symmetries. Famous examples departing from this old scheme are - apart form the high 
temperature superconductors - topological insulators, quantum critical regions connected to quan- 
tum critical points, and (fractional) quantum Hall effects, just to name a few (see [19, 20, 23-25] 
and references therein). As in particle physics, modern condensed matter theory is concerned with 
low-energy excitations that are described most efficiently by quantum field theories which reveal the 
universality of very different microscopic quantum many-body systems. However, in the strongly 
interacting cases, the "mapping" between the relevant degrees of freedom in the low-energy regime 
and the microscopic degrees of freedom are far from being clear and understood. Furthermore, it 
seems that quantum field theory alone is not enough to tackle strongly correlated systems and to 
explain these new states of (quantum) matter. Interestingly, the universality of Homes' law seems 
to go beyond the artificial distinction between traditional and modern condensed matter physics 
since it displays a relation that works for conventional superconductors and high temperature su- 
perconductors which can be regarded as representatives of the old and the yet to be developed 
framework. 

2.1 Homes' Law in Condensed Matter 

An interesting phenomenon to look for universal behavior in condensed matter systems, as adver- 
tised in the introduction to this section, is the universal scaling law for superconductors empirically 
found by Homes et al. [26] by collecting experimental results. This so-called Homes' law describes 
a relation between different quantities of conventional and unconventional superconductors, i.e. the 
superfluid density p s at zero temperature and the conductivity ctdc times the critical temperature 
T c , 

p s = Ca DC (T c )T c , (2.1) 

where Homes et al. report two different values of the constant C in units [cm] -2 for the different 
cases considered in [27]. The value C = 35 is true for in-plane cuprates and elemental BCS super- 
conductors, whereas for the cuprates along the c-axis and the dirty limit BCS superconductors they 
find C = 65. The superfluid density p s is a measure for the number of particles contributing to the 
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Figure 1: Schematic plots of the optical conductivity above the critical temperature T c and near the 
absolute zero T = for a dirty BCS superconductor. In the superconducting phase a gap develops for 
frequencies ui < 2 A. Note that in the dirty limit the quasi- particle scattering rate 1 / T is larger than 2 A. 
The missing area (red shaded region) i.e. the difference between the area under the curve of Re <t(lj, T « T c ) 
and Re a (cj, T <C T c ) which condenses into the 5-peak at ui — and thus the superfluid strength (being the 
coefficient of the J(u)-function) is proportional to the missing area. Following the definitions given in (2.6) 
we see that the area under the Kea(uj,T <C T c ) curve determines N s whereas the red shaded area yields 
N n -N s . 



superfluid phase. It can be thought of as the square of the plasma frequency 2 of the superconduct- 
ing phase Up s , because the superconductor becomes "transparent" for electromagnetic waves with 
frequencies larger than 

_ 2 4nn s e 2 
Ps = w Ps = — ■ (2.2) 

TO* 

Here n s denotes the superconducting charge carrier density which describes the number of super- 
conducting charge carriers per volume (and is very different from the superfluid density p s ), e is 
the elementary charge and to* is the effective mass of the charge carrier renormalized due to inter- 
actions. Another way to think about the superfluid density p s is the London penetration depth Al 
which is basically the inverse of the superconducting plasma frequency, such that frequencies larger 
than wp s correspond to length scales smaller than Al, i.e. p s = A^ 2 . The critical temperature T c 
is determined by the onset of superconductivity. The conductivity odc an d the superconducting 
plasma frequency wp s are for instance obtained from reflectance measurements by extrapolation to 
the oj — > limit of the complex optical conductivity cr(w), 

er D c = lim Recr(w), Up B = lim ( - cj 2 Ree(w)), (2.3) 

since the high frequency limit of the real part of the dielectric function e(uj) is given by 3 

ReeM = e 0O - a %. (2.4) 
where eoo is set by the screening due to interband transitions. 

Alternatively, the superconducting plasma frequency may be obtained from the optical conductivity 

2 The definition of the plasma frequency and its relations to the dielectric function and superconductivity is 
discussed in detail in Appendix A. 

3 A clear introduction to linear response, sum rules and the Kramers-Kronig relations 
in condensed matter can be found in [28]. 
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measured above and below the critical temperature with the help of the oscillator strength sum rule 



dcjRccr(w), (2.5) 







for the optical conductivity. This gives rise to an alternative definition of the superfluid density 
as compared to (2.3). We define the spectral weight in the normal and superconducting phase as 
follows, 



N n = / duRea(u) =^|a, N s = / 

JO T>T„ ° ^0+ 



do; Recr(w) 

t>t c ° ^0+ 



(2. 

T<T C 



The superfluid density p s describes the degrees of freedom in the superconducting phase which 
have condensed into a Dirac S-peak at zero frequency, where p s can be viewed as the coefficient 
of 8(oj). This (S-peak in the real part of the conductivity gives rise to an infinite DC conductivity 
or zero resistivity. The superfluid density is equal to the difference between the integral over the 
optical conductivity (2.5) evaluated for T < T c and T > T c and generally yields identical values as 
compared to (2.3). Using the definitions of the spectral weight (2.6) we find 

Ps = 8 (N n - N s ) . (2.7) 

This is the Ferrell- Glover- Tinkham sum rule. Note that in the definition of N s we have excluded 
the S-peak at u> = because the oscillator strength sum rule (2.5) requires that the area under 
the optical conductivity curve is identical above and below T c , i.e. in the superconducting and 
normal state. Thus (2.7) determines the missing area of the spectral weight that condensed into 
the J-peak at ui = 0, as illustrated in Figure 1 . Although the gap describes the creation of Cooper 
pairs, it is really the missing area which gives rise to superconductivity, since according to (2.7) the 
missing area is equal to the degrees of freedom which condense at zero frequency, thus forming a 
new (coherent) macroscopic ground state with off-diagonal long range order. Semiconductors for 
instance are systems exhibiting an energy gap in their spectrum as well, but are not necessarily 
superconducting since there is no missing area and hence no new ground state let alone a phase 
transition. On the other hand, superconductors retain their properties even if the energy gap is 
removed (e.g. by magnetic impurities) due to the missing area under the Re cr(uS) curve. As a caveat, 
let us note that high temperature superconductors may not satisfy the Ferrell-Glover-Tinkham sum 
rule, while it is expected to hold for dirty BCS superconductors 4 . 



In order to see clearly the relation between superfluid density and the product of the conductivity 
at the critical temperature and the critical temperature expressed by Homes' law in (2.1), we have 
reproduced Table I from [27] with and without the elemental superconductors niobium Nb and lead 
Pb in Figure 2. Furthermore, in [21,26,27,29] some possible explanations are given concerning the 
origin of Homes' law: Conventional dirty-limit superconductors, marginal Fermi-liquid behavior, for 
cuprates a Josephson coupling along the c-axis or unitary-limit impurity scattering. Furthermore, 
the authors discuss limits where the relation breaks down, which is true in the overdoped region 
of cuprates. For dirty limit BCS superconductors, Homes' law can be explained by the very broad 
Drude-peak 5 which is condensing into the superconducting <S-peak at uj — 0. The spectral weight 
of the condensate may then be estimated by an approximate rectangle of area p s w erne ' 2A in an 

4 Experimentally it is not possible to "integrate" up to ui — ¥ oo since a measurement cannot be done at arbitrary 
high frequencies, so in reality we need to introduce a cut-off frequency u) c . For high temperature superconductors 
this cut-off frequency may be higher than the experimentally accessible frequencies and thus these superconductors 
may not satisfy the Ferrell-Glover-Tinkham sum rule (see also [27]). 

5 The Drude peak is located at zero frequency where the real part of <j(uj) reaches its global maximum, see for 
example Figure f (and for more details on the Drude model see Appendix B). 
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Figure 2: Plots of the complete data in Table 1 in [27]. The error bars are calculated by A(wp) = 2a»pAo;p 
and A(oncTc) = T c Aauc- The left plot shows data from high T c superconductors and for Ba 1 _ x K x Bi0 3 , 
while the one on the right includes three data points from elemental superconductors. As shown on the 
right the elemental conventional superconductors, two data points for Nb and one coming from the Pb 
superconductor, actually give Homes' law as stated in [27]. If these three data points are ignored, as shown 
in the left panel, the linear fit is shifted and the data points are below the Homes' law line. 



optical conductivity plot, similar to Figure 1, where the gap in the energy of the superconducting 
state is denoted by 2A and ctdc is the maximum of the curve at uj = 0. According to the BCS 
model, the energy gap in the superconducting phase is proportional to the critical temperature T c 
and thus p s oc <Jt>cT c - For high T c temperature superconductors the most striking argument can 
be found in [21] which links the universal behavior to the "Planckian dissipation" giving rise to a 
perfect fluid description of the "strange metal phase" with possible universal behavior, comparable 
to the viscosity of the quark-gluon plasma. The argument, reproduced here for completeness, relies 
on the fact that the right structure of p s , ctdc and T c may be worked out by dimensional analysis: 
First, as already stated above (2.2) the superfluid density must be proportional to the density of the 
charge carriers in the superconducting state. The natural dimension for this quantity is (time)" 2 
so the product of ctdc an d T c should have the same physical dimension. Second, the normal state 
possesses two relevant time scales, the normal state plasma frequency ujp n and the relaxation time 
scale t, which describes the dissipation of internal energy into entropy by inelastic scattering. One 
of the simplest combinations is the product of the two time scales which will yield the dimension 
(time) -1 . Therefore we may take the optical conductivity to be of the Drude-Sommerfeld form (see 
Appendix B for details) given by 



"Pn' 
47T 



771 



(2.8) 



The last and most crucial step is to convert the critical temperature into the dimension (time) . 
Energy and time are related by Heisenberg's uncertainty principle and thus quantum physics and 
the idea of "Planckian dissipation" will enter, 



T h (T c y 



h 



(2.9) 
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This time scale is the lowest possible dissipation time for a given temperature. For smaller time 
scales the system will only allow for quantum mechanical dissipationless motion. Interestingly, 
at finite temperature the lower bound can only be reached if the system is in a quantum critical 
state [19]. This implies that high T c superconductors exhibit a quantum critical region above 
the superconducting dome, which is supported by experimental evidence [20,30]. To connect the 
expression on the right-hand-side of (2.1) with the left-hand-side, we can employ Tanner's law [31] 
which states that 

n s «-n n , (2.10) 

relating the superfluid density to the normal state plasma frequency, as can be seen from (2.2) and 
(2.8). This "explanation" of Homes' law will guide us in order to find a holographic realization. We 
will expand on this idea in Section 2.2. 

2.2 Homes' Law in Holography 

An obstacle towards checking Homes' law directly within holography is that due to the conformal 
symmetry and the absence of a lattice, the Drude peak of the conductivity is given by a delta 
distribution even at finite temperatures above the critical temperature, i.e. 

Recr(w) - 8{uS) Ima(cj)~-. (2.11) 

Thus, it is not possible to evaluate p s directly, which is related only to the superconducting degrees 
of freedom condensing at u> = 0. We therefore rewrite Homes' law in such a way that it becomes 
accessible to simple models of holography. As we now discuss, this can be achieved by using the 
idea of Planckian dissipation following [21] for high temperature superconductors, as outlined above 
at the end of Section 2.1, or by employing the Ferrell-Glover-Tinkham sum rule (2.7) for dirty BCS 
superconductors as can be found in [29]. 



For simplicity, let us make two assumptions about holographic superconductors: First, let us assume 
that they satisfy the sum rule which requires that the area under the optical conductivity curve is 
identical in the superconducting phase and the normal phase. Second, we assume that all degrees 
of freedom condense in the superconducting state. Using the definition of the spectral weight in 
the normal phase N n (2.6) and the definition of the superconducting plasma frequency (2.2), we 
see that both plasma frequencies must be equal 

which implies that N s = in (2.7). Here we clearly neglect possible missing spectral weight as 
explained in the second paragraph of Section 2.1. In the holographic context, similar sum rules have 
been investigated in [32-34]. Alternatively, we may assume that the holographic superconductors 
obey Tanner's law (2.13) 

n s = Bn n , (2.13) 

where n s and n n denotes the charge carrier density in the superconducting and the normal phase, 
respectively, and B is a numerical constant. With either one of these two assumptions, we may 
rewrite Homes' law in such a way that it becomes accessible to holography. Let us begin with the 
assumption that holographic superconductors fulfill the sum rule and that all degrees of freedom 
are participating in the superconducting phase: Starting from 

p s = ^ s = Ca DC (T c )T c , (2.14) 
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and assuming that the Drude-Sommerfeld model is still a useful approximation (see Appendix B 
for details) we can replace the conductivity by a characteristic time scale and the plasma frequency 
in the normal conducting phase 



and thus (2.14) reads 



v DC = (2.15) 



where r c denotes t(T c ). (2.16) can be simplified using the assumption (2.12), 

47T 

T ^ = ~C- ( 2 - 17 ) 

(2.17) holds if there is no missing spectral weight and that all the spectral weight associated with 
the charge carriers condenses in the superconducting phase and hence contributes to the 5-peak. 
This enables us to identify the plasma frequencies in the two different phases. Alternatively, we can 
use the assumption that the charge carrier densities in both states are proportional to each other. 
Starting again with 

2 

TL 6 

p B = 47r^V = Ca DC {T c )T c , (2.18) 
and inserting the Drude-Sommerfeld optical conductivity in terms of the charge carrier density n n 

<tdc = ^-r, (2.19) 

m* 

we may write (2.18) as 

C 

n s = —n n T c T c . (2.20) 

47T 

Assuming that the holographic superconductors obey (2.13) we arrive at 

r c T c = (2.21) 

Note that the proportionality constants in (2.17) and in (2.21) may not coincide. If they do not, this 
will indicate that holographic superconductors behave either more like in-plane high temperature 
superconductors or like dirty-limit BCS superconductors. In any case the above simplifications 
allow us to circumvent the need to calculate spectral weights in the superconducting phase or the 
plasma frequency in either phase (some obstructions are discussed in Section 2.3), and to perform 
the calculation solely in the normal phase. Therefore we will extract the time scales in the normal 
phase of the s- and p-wave superconductors from diffusion constants, basically the momentum and 
R-charge diffusion denoted by Dm and -Dr, respectively. In particular, since D(T C ) oc r c we obtain 

D(T C )T C = const. (2.22) 

This relation is directly accessible to holography. In fact, without including the backreaction of the 
gauge and matter fields on gravity, the diffusion constants are given by [1,35,36] 

Dm = 4^T' D *=^TdU> (2 ' 23) 
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such that 

D M T = — = const., D R T = J__ C L_ = con st., (2.24) 

47r 47r a — 2 

where d denotes the dimensionality of the spacetime, and thus Homes' law is trivially satisfied in 
this case. This can be derived simply by dimensional analysis as is done in [37]. Extending the 
calculation of the diffusion constants to include backreaction we checked analytically and numeri- 
cally that our results reduce to the known results in the limit where the backreaction vanishes, (see 
Sections 3.5 and 4.3 for details). 

2.3 The Drude-Sommerfeld Model and Holography 

The Drude-Sommerfeld model describes the properties of metals (i.a. electrical/thermal conduc- 
tivities, heat capacities) by a simple model incorporating the behavior of the underlying charge 
carriers. The optical conductivity of the Drude-Sommerfeld model (2.8) 

auc = —. — = — • (2.25) 

47T TYl 

possesses two scales, the plasma frequency Wp n setting the scale above which electromagnetic waves 
can enter the metal and the relaxation time scale t, connected to the mean free path I = v-pT of 
the charge carriers (where the charge carrier density is denoted by n). Interactions between the 
charge carriers are included by the renormalization procedure generating effective masses m* and 
correction to the relaxation time. The inverse of the real part of the optical conductivity describes 
absorption processes/loss of energy which is related to dissipation, whereas the imaginary part 
describes dispersive processes/change of phase which is related to energy transport. The resistivity 
is defined as (Rc<t(lu = 0)) _1 and is the inverse of the maximum of the real part, the so called 
Drude peak. The width of the Drude peak is set by the relaxation rate t . The imaginary part of 
(2.25) is a Lorentzian-shaped curve with maximal energy transport at resonance for lu = t . For 
very high frequencies u > wp„ the charge carriers are not able to follow the external excitation and 
thus the optical conductivity must approach zero. 



Comparing our holographic setup to the above condensed matter discussion, we first notice that we 
do not have an underlying lattice. This implies that the Drude peak turns into a delta distribution 
(or (5-peak) at uj = since momentum conservation does not allow any dissipative process. In 
the large frequency limit, the optical conductivity approaches the non- vanishing conformal result, 
i.e. a(uj — > oo) = const., since the conformal symmetry does not permit any scale. In this sense there 
is no well-defined plasma frequency above which the system will become transparent, but a finite 
absorption due to the non- vanishing real part of the optical conductivity will be retained. A possible 
workaround is to define a regulated plasma frequency [32] which obey the Kramers-Kronig relations 
and thus the sum-rules. Moreover, we cannot compute effective masses m* or more generally, we 
cannot describe any physical parameter related to the lattice. Additionally, interesting physics is 
hidden in the zero frequency limit of the optical conductivity, which cannot be resolved without a 
lattice due to the (5-peak arising from momentum conservation 6 . One way to remedy these problems 
is to introduce a lattice explicitly as is done in different ways in [17,18,38,39]. One of the results 
in [17] is the emergence of the Drude-Sommerfeld conductivity (B.6) in the low-frequency limit and 
- even more excitingly - in the mid-frequency range they find scaling behavior similar to cuprate 
superconductors. The use of an explicit lattice may be avoided by using the simplified form of 

6 Technically, momentum conservation can be broken by neglecting the effects of backreaction onto the geometry 
on the gravity side. Physically this is not very helpful since the fixed geometry introduces an artificial dissipative 
reservoir. 
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Homes' law stated in (2.17) or (2.21) which assumes either the validity of the sum rule or Tanner's 
law, respectively. This fits in our overall scheme of looking at universal scaling behavior which 
should not be affected by an underlying microscopic lattice. 



3. Holographic s-Wave Superconductor 

In this section we give a self-contained exposition of the solution to the equations of motion derived 
from the Einstein-Maxwell action minimally coupled to a charged scalar field. For obtaining a family 
of holographic s-wave superconductors, we include the backreaction of the gauge field and the scalar 
field on the metric into our analysis. After explaining the normal phase solutions, i.e. solutions with 
a vanishing scalar field, we will employ a quasi-normal-mode analysis in order to determine the onset 
of the spontaneous condensation of the operator dual to the charged scalar field. For this purpose, 
we expand the action up to the quadratic order in the fluctuations about the background and 
derive the corresponding equations of motion. Finally, we determine the R-charge and momentum 
diffusion related to the gauge field fluctuations and the metric fluctuations, respectively. This allows 
us to calculate the function relevant for testing Homes' law at finite backreaction. 

3.1 Einstein-Maxwell Action 

The best known (bottom up) holographic model to describe superconductivity is given by the 
Einstein-Maxwell action coupled to a charged scalar on the gravity side dual to a field theory with 
a conserved £7(1) current and an operator describing the condensate [13,40,41]. Rescaling the gauge 
field A — > eA and the scalar field $ — > e* allows us to write the action in the form 



2k 2 

where we have a dimensionless coupling parameter 



R-2A-^ r (^F ab F ab - |V* - U$| 2 - V(\$\) 



(3.1) 



2 



a'L z = ^, (3.2) 

describing the strength of the backreaction onto the geometry exerted by the gauge field and the 
scalar field. On the field theory side, this parameter describes the ratio of charged degrees of freedom 
to the total degrees of freedom and thus can be considered as an effective chemical potential or in a 
loose sense some kind of "doping". The equations of motions for the charged scalar field <j> derived 
from (3.1) are given by 



(V„ + iA a ) (V a + iA a ) <r - = o, 

(V -iA a )(V°-iA a )$- ±v'Q$\y* 



(3.3) 



2 " "|*| 

Additionally, we have the Maxwell equations in curved spacetime 

V a F ab = j b , (3.4) 

which can be simplified to 

-J=9d,d c (V=9F cb ) = i[$* (V 6 - iA b ) $ - $ (V 6 + xA b ) <T ] . (3.5) 
Finally, the Einstein equations sourced by the gauge field and by the scalar field are 

Rab - \Rgab + A 3ab = a 2 L 2 (IX + l£) . (3-6) 
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with 

Tab 1 = 9 cd FacFbd - ^9abF c dF cd , (3.7) 

T? b = (V Q $* + L4 a $*) (V h $ - L4 h $) + (V a $ - L4 a $) (V & $* + iA b $*) 

- 5ab (V c $* + L4 C $*) (V c $ - iA c $) - g ab V(\$\). (3.8) 

The most general planar and rotational symmetric solution to the set of equations (3.3)-(3.6) is the 
AdS-Reissner-Nordstrom black hole solution with scalar hair . Note that it is sufficient to assume 
that only quadratic terms are present in the potential V(|$|) = m 2 l^ 2 , since we are only interested 
in the behavior near the critical point where higher order interactions do not contribute. 

3.2 Normal Phase Solutions of the Background Fields 

As we are only interested in the phase diagram and since we are including the backreaction. it is 
advisable to simplify our discussion by considering the normal phase background solutions, i.e. $ = 
0, which reduces the equations of motion to 

Rab - \Rgab ~ ^^9ab = O? L?T™ , (3.9) 

g ab V a F bc = 0. (3.10) 

The solutions to these equations are given by 



ds 2 = ^(-f(u)dt 2 + dx 2 +^) 

u 1 v .f{u); 



with 



f{v) = l-M-^ + Q'-^-^, M = l + Q\ Q= ^— fIUH ^ (3.12) 



*H M H 

where M and Q denotes the dimensionless mass and the dimensionless charge of the black hole, 
respectively. The gauge field solution including the constraints on the black hole horizon may be 
written as 

A = fi^l-^jdt, (3.13) 

where /i is the chemical potential, describing the difference in the electric potential between the 
horizon of the black hole at u = i*h and the boundary of the AdS space (at u = 0) and thus 
the energy of adding another charged particle to the black hole. Furthermore we may define a 
temperature related to the Bekenstein entropy by 



2tt 



1 d 

V-9u 



fg~^L du 



(3.14) 

47Tli H 



Since we are dealing with a scale invariant theory in the UV, the only dimensionless physical 
parameter is the ratio T /pL which will be controlled (numerically) by the dimensionless chemical 
potential p, — /iuh- Using (3.12) we may write this ratio in terms of the charge of the black hole, 

T _ _ t-v-w 

fl 4:77/1 4tt _ I d-l Q 
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d = 3 



d = 4 
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a a 



Figure 3: Phase diagram of the holographic s-wave superconductor for d — 4 and d = 3 as a function of the 
backreaction parameter a, depending on the scalar field mass. Colored regions (for d = 3 are encoded as 
m 2 L 2 = 4, m 2 L 2 = 0, m 2 L 2 = -2 and for d = 4 we have m 2 L 2 = 5, m 2 L 2 = 0, m 2 L 2 = -3, m 2 L 2 = -4) 
show phases where the scalar field condenses yielding a superfiuid phase. 



and conversely 



Q± = 



-2it^d-l (£) ± W47T 2 (rf- 1) 



d(d-2) 2 



a(d-2)i 



(3.16) 



In order to be consistent with the probe limit a — > 0, we pick the positive solution of the quadratic 
equation for the black hole charge Q, such that limo^o Q+ = 0. 

3.3 Quasi-Normal-Mode Analysis and Phase Diagram 

We use the holographic dissipation-fluctuation theorem [42] in order to determine the complex- 
valued Green's function of scalar fluctuations about the fixed scalar background in the normal 
phase. Therefore we derive the quadratic action in the fluctuations 



</>(£, x, u) = + 54>(t,x, u), 
A a (t,x,u) = A*(u) + a a (t,x,u), 
g a b(t,x,u) = G ab {u) + h a b(t,x,u). 



(3.17) 



The detailed calculations and solutions are shown in Appendix C. For convenience we state the 
equation of motion for the scalar field fluctuations (C.5), 



b"{u) + 



/(«) 



i 



54>'{u) 



f(u) 2 



L 2 m 2 



f(u) u 2 f(u) 



5cj){u) = 0, 



(3.18) 



which we will use in order to determine the phase diagram. We are interested in the spontaneous 
symmetry breaking induced by the condensation of the operator dual to the scalar field <f>. For a 
second order phase transition we can use the dissipation-fluctuation theorem to look for instabilities 
in the fluctuations. Thus, for a given fixed value of the backreaction a, we vary the dimensionless 
parameter T /pL and look for a quasi-normal mode crossing the origin of the complex frequency 
plane. Numerically, we solve (3.18) with infalling boundary conditions at the black hole horizon 
and calculate the corresponding retarded Green's function for k = and ui = 0. We expect to 
find a gapless mode which is related to a pole in the retarded Green's function, such that we may 
determine the set of parameters (a, t /m) yielding critical curves as shown in Figure 3. Additionally, 
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we can also vary the mass of the scalar field related to the scaling dimension of the dual operator 
via 



A± = i (d± \/d 2 + 4L 2 m 2 ^j , m 2 L 2 = A± (A± - d) . (3.19) 

The critical value of the backreaction for T = can be determined analytically by looking at the 
IR behavior near the black hole horizon. The corresponding charge of the black hole is given by 

(3.20) 



d-2 



Inserting this value into the metric (3.11), we find a double zero giving rise to an AdS2 x H d 1 
metric and thus an IR fixed point 

d4 dS2 =L 2 (Vdi+^), (3.21) 

with g = (u — uh) and rescaled t and x coordinates, accordingly. The AdS2 radius is related to the 
AdSd+i radius by 

2 L 2 

LAdS ? ~ djd^Ty (3 - 22) 

The near horizon expansion of the background gauge field reads 

2 djd - 1) 2 d(d-l) 2 

Finally, (3.18) can be written as 

S^(g) + %'(g) - 6tP(q) = 0, (3.24) 

and the effective mass term can be compared to the Breitenlohner-Freedman stability bound for 
AdS 2 , 

2 2 1 I t2 2 M X r2 ° 



L AdS 2 m eS = d t d _ l \ [ L m ~o^)- L AdS 2 m BF = ( 3 - 2 5) 



d(d-l) V a 2 
This leads to a condition on a for scalar field condensation to occur, 

1 



L 2 



i))- 



We see if the mass is already below the AdS2 Breitenlohner-Freedman bound, there is no critical 
value for a and hence no quantum critical point or phase transition at zero temperature between 
the condensed phase and the normal phase. The different masses used in the numerical calculation 
and the corresponding values for the scaling and the critical backreaction strength are listed in 
Table 1. Our results as displayed in Figure 3 show that the critical temperature decreases with 
increasing backreaction strength a. Moreover, if the scalar mass is larger than a critical value, the 
critical temperature goes to zero for a finite value of a. This is the case most reminiscent of a real 
high T c superconductor, when the dome in Figure 3 has similarities with the right hand side of the 
dome in the phase diagram of a high T c superconductor. 



The physical interpretation of a is that it corresponds to the ratio of the number of SU(2) charged 
degrees of freedom over all degrees of freedom [43]. The phase diagrams above indicate that an 
increase of a reduces the numbers of degrees of freedom which can participate in pair formation 
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Table 1: List of the critical value for the strength of the backreaction a for different masses in three and four 
dimensions. Note that the instability condition is satisfied for a < a c . In particular for the negative values 
we do not have a critical value of a 6 E, so in this case for T = we always find a condensed/superfluid 
phase. 

and condensation, such that T c is lowered. A similar mechanism also seems to be at work when 
adding a double trace deformation to the holographic superconductor [44] , and has been discussed 
within condensed matter physics using a BCS approach in [45]. For holographic superconductors, 
this mechanism is most clearly visible for the top-down holographic superconductors involving D7 
brane probes [46,47] in which the dual field theory Lagrangian and thus the field content of the 
condensing operator are known. In these models, there is an U (2) symmetry which factorizes into 
an SU(2)i x U(1)b, i-e. into an isospin and a baryonic symmetry. A chemical potential is switched 
on for the SU(2)i isospin symmetry and the condensate is of /o-meson form, 

Jx 1 = i>a l lxip + §v X dx<P = $fj x ipi + i>i7 x ^ + bosons, (3.27) 

where ip = (ip-\, ip±) and <p = <p\J) are the quark and squark doublets, respectively, which involve 
up and down flavors, with a 1 the Pauli matrices in isospin space and 7^ the Dirac matrices. As an 
additional control parameter, a chemical potential ps for the baryonic U(1)b symmetry may be 
turned on. This leads to a decrease of T c [48] (see also [49]) which may be understood as follows: 
Under the U(1)b symmetry, ip and ip have opposite charge. The same applies to <p and (p. Turning 
on hb leads to an excess of ip over ip degrees of freedom, which implies that less degrees of freedom 
are available to form the pairs (3.27). The same applies to (p and <p degrees of freedom. Thus in 
this case, charge carriers in the normal state are also present in the superconducting phase, leading 
to the formation of a pseudo-gap. 

3.4 Diffusion Constants 

As explained in Section 2, an analysis of Homes' law requires the definition of a relevant time scale. 
Possible candidates for time scales are associated with diffusive processes. We may determine two 
different types of diffusion, momentum diffusion Dm and R-charge diffusion Dr, respectively. The 
former can be related to the shear viscosity by studying hydrodynamic modes, 

D M = -Jp, (3.28) 

such that it can be solely described by thermodynamic quantities using the famous result (1.1). 
The latter, however, cannot be calculated solely by thermodynamic quantities in general, such that 
an independent calculation for each background is necessary. Usually the charge diffusion constant 
can be read off from the dispersion relation w = — iDk 2 of the conserved current, which can be 
derived from Fick's law and the continuity equation 

j = -DVp, and p + V-j = 0, — > p-D\7 2 p = Q. (3.29) 

Fortunately, this work has already been done by [35] without backreaction. 
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Momentum Diffusion with Backreaction 

Since we are working with a fixed chemical potential, we take the grand canonical ensemble to 
describe the momentum diffusion in the thermodynamic limit. Using the gauge/gravity duality the 
grand canonical potential is given by 

= = TI, (3.30) 
where / denotes the regularized Euclidean on-shell action given by 

7 =-!?¥4 (i+<52) ' 

Using the Gauge/Gravity dictionary, we derive the charge density, the energy density and the 
entropy density from (3.30), with the help of the relations 



4tt L 



d-l 



2k 2 ut' 1 



-T^T^ (3-32) 
n = ^{d-l)(d-2) Ld ~ 1 a ® 



^ ut 1 ' 



The black hole horizon «h should be understood as a function of /i or T, respectively, defined by 
solving either (3.12) or (3.14) for un, whereas Q is a function of T /n defined by the inversion of 
(3.16). Starting from the definition of the momentum diffusion and using the relations (1.1) and 
the first law of thermodynamics e = Ts — P + /in, we find 



v s 1 1 

P Ts + ^n AttTI , 



D M - 7~~i 5 - — - JZ^ , m n (3-33) 



Inserting the relations (3.15) for m/t and (3.32) for n /s we finally arrive at 



D u -*(l + «<- 1 *' ) . (3.34) 



AttT \ d-{d-2)Q 2 
R-Charge Diffusion with Backreaction 

As explained in the previous subsection it is not straightforward to generalize the R-charge diffusion 
with backreaction to arbitrary dimensions. To our knowledge the functional dependence on the 
backreaction of in d = 3 dimensions is not known yet, let alone for arbitrary dimensions d. 
Surely, this would be an interesting result which is beyond the scope of this work. In d = 4 the 
R-charge diffusion with backreaction is determined in [50]. Converting the expression for Dr into 
the form given by our conventions we find 

_ 1 (2-Q2)(2 + Q2) 
DR -±nf 2(1 + ' ( 3 - 35 ) 

where the dimensionless black hole charge is defined in (3.12). Furthermore, we may look at the 
ratio of the two different diffusion constants, 

^=2 + Q 2 . (3.36) 
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Figure 4: Plots of momentum diffusion constant Dm for d = 3 (left panel) and R-charge diffusion Dr 
and momentum diffusion Dm for d — 4 (right panel) versus the backreaction strength a. The color coding 
of the different scalar field masses follows Figure 3. The probe limit properties for a — are reproduced, 
i.e. for d — 3 we find Dm = 1 /a and in the case of d — 4 Dm = V 4 an d Dr = 1 /2, independently of the scalar 
field mass. 



3.5 Discussion 

Diffusion without backreaction for arbitrary d 

First of all we check our results by taking the limit of vanishing backreaction. i.e. setting a and 
Q to zero in (3.34) and (3.35). We see that the dimensionless temperature T (3.14) has the fixed 
value of d /i-n and thus 

D M = = _ (3.37) 
is equal for all masses of the scalar field. Similarly, the R-charge diffusion is given by [37] 

D «=^Td^ = d^2> 



and the ratio is given by 



D M d-2' 



(3.39) 



Comparing these results to our numerical solutions shown in Figure 4, we see that for a = 0, we 
obtain for d = 3 dimensions a dimensionless value of 1 /3 for Dm- This is consistent with the values 
obtained form the analytic calculation without backreaction (3.37). In particular, we see that the 
momentum diffusion does not depend on the mass of the background scalar field since all diffusion 
constants for different masses (indicated in the figure by different colors) converge to the same value. 
Similarly for d — 4, we check our numerical values for the momentum and the R-charge diffusion by 
comparing to the results without backreaction. As is displayed in Figure 4, we find that Dm = l /i 
and Z>r = Y 2 as stated in (3.37) and (3.38), as well as the correct value of the ratio d r/d u — 2. 
Again our numerical results are consistent with the analytic solutions without backreaction and we 
see the same convergence effect for the different masses of the scalar field for each diffusion constant 
displaying the independence on the scalar fields' mass. 

Momentum Diffusion for d = 3 

Inserting the critical values for at a fixed a, we may check how much the momentum diffusion 
is varying with respect to a. As shown in Figure 5, there is a mild variation of the momentum 
diffusion in the sense that we are approaching the values for a — in a linear way. In order to 
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Figure 5: The constant Cm related to momentum diffusion Dm for d — 3 and different masses of the 
scalar field (color coded as m 2 L 2 = 4, m 2 L 2 = and m 2 L 2 = — 2) plotted depending on the strength of 
the backreaction a, defined in (3.2). The left panel shows the numerical values found by a minimization 
algorithm up to the a where poles of higher excitations in the Green's function overlays with the lowest 
modes, thus leading to a breakdown of the algorithm. Only in the case of m 2 L 2 = —2, we have tested our 
algorithm up to a = 1. In the right panel we included the analytically determined critical points a c for 
in L 2 = 4 and m 2 L 2 = given in Table 1 which corresponds to the zero temperature case with Q 2 — d /(d-2) 
and thus Cm = 0. Due to the reduced number of data points we use an interpolation function in the region 
a 6 [0.35, sf^fi ~ 0.82]. For this reason the curves for m 2 L 2 = 4 and m 2 L 2 = should not be trusted to be 
very accurate. Since for m 2 L 2 — —2 there is no critical value of a, so Cm approaches zero asymptotically. 



compare our numerical results to the constant governing Homes' law rewritten in the form (2.17) 
or (2.21), we introduce the function 

C M = ^T C D M (T C ) = (l + 7~^j , (3-40) 

plotted in Figure 5. Here we assume that the diffusive process is proportional to a time scale 
t c oc Dm(T c ), such that the proportionality does not depend on Q or a. 

R-Charge Diffusion and Momentum Diffusion for d = 4 

We now repeat the discussion of the previous subsection for the four-dimensional case. In addition 
to the momentum diffusion, we may also calculate the R-charge diffusion which will give us an 
additional time scale to choose. The constant related to 4wt c T c for the R-charge diffusion reads 

_ (2-Q 2 )(2 + Q 2 ) 
Cr ~ 2(1 + Q») ' (3 ' 41) 

Note that this function is running from 2 to for a € [0, a c ] since the R-charge diffusion without 
backreaction is given by 1 /(2kT) see (3.38). For d = 4, the constant for the momentum diffusion 
reads 

CM=(l + ^)"\ (3.42) 

Note also that the ratio c ^/c m is given by 2 in compliance with the ratio for the diffusion constants 
(3.39). The plots of Cm and Cr are shown in Figure 6 and Figure 7, respectively. We see that the 
numerically solutions are virtually the same as in the three dimensional case. In both dimensions 
the numerical solutions display the same properties: 
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Figure 6: For d = 4 dimensions: Dimensionless momentum diffusion given by Cm depending on the 
backreaction a (for different masses of the scalar field color coded as m 2 L 2 — 5, m 2 L 2 = 0, m 2 L 2 — —3 and 
m 2 L 2 = —4). As already explained in Figure 5, the endpoint in the curves for scalar field mass rn 2 L 2 — 5 
and m 2 L 2 = as shown in the right panel is set by the critical values of a listed in Table 1. Due to the 
interpolation these curves should be taken with a grain of salt in the vicinity of the critical point a c . 



d=4 d=4 
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Figure 7: Dimensionless R-charge diffusion given by Cr = 47tT c Dr. Note that we have used the same 
prefactor in the definition (3.4f) as for the dimensionless momentum Cm, which leads to the probe limit 
value of two. The color coding of the scalar field masses is identical to Figure 6. 

• Instead of the behavior D(T C )T C — const, predicted by Homes' law in combination with the 
sum rules as discussed in Section 2.2, we observe a decrease of this quantity with increasing 
backreaction strength a. We discuss possible explanations for this behavior below in section 
5. 

• A test on our calculations is provided by the fact that numerically, we reproduce the re- 
sults known from probe limit calculations where by definition Cm and Cr are one or two, 
respectively. 

• The region where our minimization algorithm is working safely is approximately given by 
a w 0.35. For higher values we use an interpolation function for the curves associated to the 
non-negative scalar field masses including the zero temperature values obtained analytically. 

• At the endpoints where we encounter a quantum critical point with vanishing critical temper- 
ature, Cm and Cr are zero since these points correspond to the critical a stated in Table 1. 
For the negative masses there exists no critical a and the curves are reaching zero only asymp- 
totically. 
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4. Holographic p-Wave Superconductor 

In this section we discuss the Einstein- Yang-Mills theory including backreaction in asymptotically 
AdS§ which give rise to a holographic p-wave superfluid at low temperatures [43]. In addition to 
the Abelian symmetry the rotational symmetry 50(3) is spontaneously broken down to SO(2). As 
far as Homes' law is concerned, we observe very similar dependence of D(T C )T C on the backreaction 
as in the s-wave case. 

4.1 Einstein- Yang-Mills Action 

In addition to gravity, we consider here gauge fields transforming under some gauge group. This 
system is described by the Einstein- Yang-Mills action. In the following we specialize to (4 + 1)- 
dimcnsional asymptotically AdS space and to the gauge group SU(2) with field strength tensor 

F« v = d^Al - d v Al + e abc AlAl, (4.1) 

where e abc is the total antisymmetric tensor and e 123 = +1. In analogy with QCD, we call the 
SU(2) isospin symmetry. 



The Einstein- Yang-Mills action reads 



S = 



-fJ 



— (R-A) —F a F a ^ v 

2k% K > An 2 

Zfi 5 "YM 



pV 



(4.2) 



where k 5 is the five-dimensional gravitational constant, A = —12/R 2 is the cosmological constant, 
with R being the AdS radius and <?ym the Yang-Mills coupling. 



The Einstein and Yang-Mills equations derived from the above action are 

Rpv + ~j^9liv — K 5 (t^v — T^TpQ^^j , 

where the Yang-Mills energy-momentum tensor T^ v is 



(4.3) 



9ym . 



pa pap 



pa paap 
y^,v r ap r 



(4.4) 



A known solution of these equations of motion is the AdS Reissner-Nordstrom black hole discussed 
in 3.2 where a = ks/^ym- The gauge field is given by 

A = n (l - \) r 3 dt . (4.5) 

where r 3 = er 3 /2i with cr 3 the third Pauli matrix. We consider the diagonal representations of the 
gauge group since we may rotate the flavor coordinates until the chemical potential lies in the third 
isospin direction. 

Another solution which corresponds to the holographic p-wave solution may be obtained if we 
choose a gauge field ansatz [43,51] 

A = 0(r)r 3 dt +w(r)r 1 dx . (4.6) 
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The motivation for this ansatz is as follows: In the field theory we introduce an isospin chemical 
potential by the the boundary values of the time components of the gauge field <f>. This breaks the 
SU(2) symmetry down to a diagonal U(l) which is generated by r 3 . We denote this U(l) as U(l)s. 
In order to study the transition to the superfluid state, we allow solutions with non-zero ( J*) such 
that we include the dual gauge field A\ = w in the gauge field ansatz. Since we consider only 
isotropic and time-independent solutions in the field theory, the gauge fields exclusively depend on 
the radial coordinate r. With this ansatz the Yang-Mills energy-momentum tensor defined in (4.4) 
is diagonal. Solutions with (J*) ^ also break the spatial rotational symmetry SO(3) down to 
SO(2) 7 such that our metric ansatz will respect only 50(2). In addition, the system is invariant 
under the Z2 parity transformation Pu : x — > —x and w — > —w. Furthermore, given that the 
Yang-Mills energy-momentum tensor is diagonal, a diagonal metric is consistent, 

ds 2 = -N(r)a(r) 2 dt 2 + — ^--dr 2 +r 2 f(r)- 4 dx 2 +r 2 f(r) 2 (dy 2 +dz 2 ) , (4.7) 
7V(r) 

with N(r) = -2m{r)/r' 2 +r 2 /R 2 . 



Inserting our ansatz into the Einstein and Yang-Mills equations leads to six equations of motion for 
m(r), cr(r), /(r), <p{r ), w(r), ip(r) and one constraint equation from the rr component of the Einstein 
equations. The dynamical equations may be written as 



a 2 rf A w 2 <f 2 r 3 a 2 <j>' 2 fr 3 f' 2 



6Na 2 6cr 2 
a 2 f 4 w 2 <j) 2 (2rf 2 a 2 fw 12 



ZrN 2 o \ f 2 3r 

a 2 f 5 w 2 (b 2 a 2 f b w' 2 ,, /3 /' N' a 



Sr 2 N 2 a 2 1 3r 2 f [ r f + N + a ) ' 

V = 2 AT ~ V 

web 2 ,/l 4/' N' a' 
N z a z \r j N a 

The equations of motion are invariant under five scaling transformations (invariant quantities are 
not shown), 

(I) a ->• Act, </>->• \<f>, 

(II) / A/, w \~ 2 w, 

4 (4-9) 

(III) r — > Xr, m — > A m, w — > Xw, <fi — > Xcj), 

(IV) r -> Xr, m X 2 m, R -> XR, (j) -> X^ 1 ^, a -> Aa, 

where in each case A is some real positive number. We use (I) and (II) to set the boundary values 
of both a and / to one, so that the metric will be asymptotically AdS. Also we can use (III) to set 
Th to one, but we will keep it as a bookkeeping device. We use (IV) to set the AdS radius R to one. 

4.2 Phase Diagram 

Let us now discuss the phase structure of this theory. At temperature below the critical temperature 
the thermodynamically favored phase is the holographic superfluid. By varying the parameter a, 
the critical temperature changes. In addition in [43] it was found that the order of the phase 



7 Note that the finite temperature and chemical potential already break the Lorentz group down to 50(3). 



- 20 - 



0.10 
0.08 
0.06 

T 

0.04 
0.02 
0.00 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

a 

Figure 8: Phase structure of the p-wave superconductor theory: In the blue and red region the broken 
phase is the thermodynamically preferred phase, while in the white region the Reissner-Nordstrom black 
hole is the ground state. In the blue region the Reissner-Nordstrom black hole is unstable and the transition 
from the white to the blue region is second order. In the red region the Reissner-Nordstrom black hole is 
still stable. The transition form the white to the red region is first order. The black dot determines the 
critical point where the order of the phase transition changes. In the green region we cannot trust our 
numerics. 

transition depends on the ratio of the coupling constants a. For a < a c — 0.365, the phase 
transition is second order while for larger values of a the transition becomes first order. The critical 
temperature decreases as we increase the parameter a. The quantitative dependence of the critical 
temperature on the parameter a is given in Figure 8. The broken phase is thermodynamically 
preferred in the blue and red region while in the white region the Reissner-Nordstrom black hole is 
favored. The Reissner-Nordstrom black hole is unstable in the blue region and the phase transition 
from the white to the blue region is second order. In the red region, the Reissner-Nordstrom black 
hole is still stable however the state with non-zero condensate is preferred. The transition from the 
white to the red region is first order. In the green region we cannot trust our numerics. At zero 
temperature, the data is obtained as described in [52,53]. 

4.3 Diffusion Constants 

We now determine the diffusion constants at the critical temperature. Since the system is still 
described by a Reissner-Nordstrom black hole as in the holographic s-wave superfiuid and thus the 
equations of motion for the fluctuations coincide, we can use the same expressions to calculate the 
diffusion constants for the holographic p-wave superfiuid. The only change is the dependence of the 
critical temperature on the parameter a. The numerical results are shown in Figure 9. Comparing 
the phase diagram of Figure 8 and the diffusion constants of Figure 9, we see that they are virtually 
identical. Furthermore, the analytical expression converting the curve of critical values of T /V into 
the product D(T C )T C is the same for both holographic superconductors. Thus, it is not surprising 
that we will find answers that have a very large resemblance. Moreover, the same is true for the 
s-wave superconductors and the comparison between the s-wave and the p-wave superconductor: 
All curves for D(T C )T C are very similar in both cases. In the following section we will give a detailed 
conclusion and explain some possible mechanisms which may lead to the decrease of D(T C )T C in the 
backreacted case. 

5. Conclusions 

Let us summarize the main results of our paper: Motivated by Homes' law we have analyzed the 
diffusion constants in holographic s- and p-wave superconductors with backreaction. In particular 
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Figure 9: Cr = AtiT c Dr(T c ) and Cm = A-kT c Dm{T c ) related to charge diffusion Dr and momentum 
diffusion Dm depending on the strength of the backreaction a for the holographic p-wave superfluid with 
d — 4. The blue line corresponds to the constants evaluated at the temperature at which the Reissner- 
Nordstrom black hole becomes unstable (see blue line in Figure 8). The red dots corresponds to the constants 
evaluated at the critical temperature at which the superfluid phase is thermodynamically preferred (see red 
dots in Figure 8). 



we have discussed the temperature and density dependence of the momentum and R-charge dif- 
fusion for various masses of the scalar field (in the s-wave case) and for different strengths of the 
backreaction. We have found that the diffusion constants decrease with increasing backreaction 
and that the decay of the diffusion constants increases with increasing mass of the scalar field. For 
non-negative masses of the s-wave superconductor scalar field, we find an emerging AdS2 at zero 
temperature which defines a critical strength of the backreaction a c . Above this critical backreac- 
tion there is no phase transition to the superfluid/condensed phase, thus implying that a c defines 
a quantum phase transition between the normal and the condensed phase. These results are of 
intrinsic interest within holography. 

Let us now discuss these results in relation to Homes' law: We have found that without backreac- 
tion, holographic superconductors obey Homes' law since t c T c = const, in all spacetime dimensions 
greater than two. In the case of s-wave superconductors, this holds for all scalar field masses. With- 
out backreaction, the result t c T c = const, is almost automatic since the diffusion constant D scales 
as 1/T. Turning on the backreaction, we find that D(T C )T C is no longer constant, but decreases 
with increasing backreaction strength. Both for s-wave and p-wave holographic superconductors, we 
obtain virtually the same curves for the various diffusion constants in three and four dimensions. In 
particular, the s-wave superconductor shows this behavior in three and four dimensions for different 
masses of the scalar field. This is a strong indication that there is a universal principle at work, 
although we do not find a true constant. There are several possibilities how corrections may arise: 



1. The simplest explanation is the fact that we cannot assume the validity of the assumption 
(2.12) for holographic superconductors. As discussed at the end of Section 3.5, increasing 
the backreaction a leads to the formation of a pseudo-gap and there are normal state charge 
carriers present in the superconducting phase as well. The dependence of the pseudo-gap on 
the backreaction has been studied holographically in [54] where it was found that the pseudo- 
gap arising in the superconducting phase becomes larger with increasing backreaction, as 
clearly visible in Figure 4 of [54]. Assuming that the Ferrell-Glover-Tinkham sum rule (2.7) 
is still valid in the presence of the pseudo-gap, N s in (2.7) is no longer zero. This leads to 
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correction terms to (2.12) and (2.17), 

c4 s =o4 n -8iV s o Tc T c = ^(l-^. (5.1) 

Thus, we would expect the "constant" in (2.21) to decrease and this is exactly what we see in 
Figure 5, 6, 7 and 9. 

2. Even more dramatically, the sum rule (2.7) may not be valid for holographic superconductors 
in the presence of the backreaction. In this case Homes' law will be implied by Tanner's law 
(2.13) which applies to high T c superconductors. It would be interesting to study Tanner's 
law in a holographic context (regardless of its relevance to Homes' law). Our result leads to 
the conjecture that the relation between the superconducting charge carrier and the normal 
state charge carrier concentration is dependent on the strength of the backreaction 

n s = B(a)n n , (5.2) 



such that (2.21) is modified by the monotonically decreasing function B(a) yielding 

r c T c = ^. (5.3) 

3. The proportionality between the time scale and the diffusion constant could in principle 
depend on a and Q, i.e. a function A(a), say. In this case the question arises if there might 
be some additional dynamics concerning diffusion in holographic superfluids/superconductors 
which would have to be taken into account. This would lead to the modified version of (2.22), 

D <™ " jgjc- < 5 ' 4 > 

4. From a condensed matter point of view, we should look at the dominant time scale which 
is given by the energy relaxation time. This needs not necessarily to be the same as the 
momentum time scale. However, the Einstein-Higgs models obey relativistic symmetries and 
therefore the momentum diffusion/time scales are related to the energy diffusion/time scales. 
In fact, the energy diffusion constant is identical to the momentum diffusion as expected 
from the symmetries, up to a factor depending on the spacetime dimensions: Conformal 
symmetry implies that the bulk viscosity vanishes, which ensures that hydrodynamic modes 
corresponding to energy and momentum diffusion are related. 8 So the natural choice is indeed 
to look at momentum and R-charge diffusive processes to determine the relevant time scales, 
as done in this paper. Nevertheless, this issue deserves further investigation. 

5. A more speculative and interesting correction could originate from the fact that we already 
have an infinite Drude peak in the normal phase indicating an ideal conductor, i.e. Re <j{uS) oc 
5 (u>), since there is no lattice present in our model. Entering the superconducting phase seems 
to add additional spectral weight due to the condensation of the scalar field. This in turn 
implies that the perfect conducting metal in the normal state does not turn into a supercon- 
ductor but that the ideal conductor and the superconductor coexist in the superconducting 
phase. Here we would again expect a violation of equal spectral weight, or of the Ferrell- 
Glover-Tinkham sum rule (2.7), as well as having a simple relation between the normal state 
charge carrier concentration and the one in the superconducting state. In order to really 



8 A natural candidate for energy diffusion is the sound attenuation V following from the dispersion relation ui(k) 
ti soun< jfc — iFk 2 . For a conformal fluid, V is determined by the momentum diffusion via V = £*M- 
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understand what happens in the superconducting phase, one needs to determine the behavior 
of p s for small temperatures and to check that only the contribution to the superconducting 
state is considered. It is believed that exactly at T = the normal state delta shaped Drude 
peak vanishes and the coefficient of the delta distribution is really the supcrfluid strength, yet 
this is not explicitly shown so far to our knowledge. 

6. Outlook 

As is explained in Section 5, the most striking explanation for corrections to t c T c = const, in the 
backreacting case comes from additional degrees of freedom in the pseudo-gap. Let us give a list of 
open questions which needs to be addressed to make further progress: 

• An explicit calculation to verify the validity of the sum rules of holographic superconductors 
is in itself an interesting question. Therefore one needs to calculate the superfiuid density 
including its dependence on the backreaction at zero temperature as well as the optical con- 
ductivity above the critical temperature and near zero temperature. The conceptual setting 
is quite clear, but the numerical implementation might be challenging. 

• In order to confirm the assumptions made in (2.13) one needs to check if the holographic 
superconductors show a relation between n s and n n , i.e. Tanner's law. The same challenges 
might arise as in the aforementioned point. 

• It would be interesting to understand if the holographic superconductors consists of an ideal 
metal coexisting with a superconductor following the speculative point 5 made in the conclu- 
sion. If so, one needs to devise a scheme to calculate the superfiuid strength by removing the 
influence of the perfectly conducting metal. There are several ways to achieve this: First of all 
it would be interesting to look for the Meifiner-Ochsenfeld effect, i.e. to calculate the transver- 
sal response of the gauge field which is only sensitive to the superfiuid strength deep in the 
superconducting phase. Alternatively, a further idea is to change the geometry such that the 
degrees of freedom of the normal state become very massive and are thus "gapped out" of the 
spectrum. Possible candidates are the hard-wall geometry, studied in the condensed matter 
context in [55] , or a more smooth realization emerging in the AdS-soliton solutions [56] , focus- 
ing in particular on the AdS-soliton to AdS-soliton superconductor transition. A suppression 
of normal low energy spectral weight has also been observed recently in a different context 
in [57]. It will also be interesting to consider holographic systems with lattice structure, such 
as [17] and [18], in the context explained here. These do not have the <5-peak in the normal 
phase, as expected for a system with an underlying lattice. 

• Since within condensed matter physics, energy diffusion is relevant to Homes' law, a system- 
atic investigation of the relation between energy diffusion and momentum diffusion in the 
holographic context appears to be highly desirable, in particular in the presence of finite 
density and backreaction. As explained in the conclusion, our assumption that energy and 
momentum diffusion are related within holography is well-motivated. Nevertheless this is a 
crucial issue which requires systematic study. 

• It would also be useful to calculate the dependence of fermionic excitations on the backreac- 
tion, for instance by generalizing the work of [53] and [58]. 

In any case it is very important to disentangle the different contributions to the 5-peak in the 
superconducting phase in order to learn the differences which arise in holographic superfluids/super- 
conductors as compared to "real world" systems. 
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Appendices 

A. Plasma-Frequency 

Here we give a brief overview over the plasma frequency and the relations to the dielectric function 
and its role in terms of the superfluid density. Formally, the plasma frequency is defined as 



/>oo 

/ dwRccr(cj), (A.l) 
Jo 



using the optical sum rule. Using the relations between the complex dielectric function e(u>) and 
the complex optical conductivity a(ui), 

47ri uj 

e = £co H cr(u;) Recrfw) = — Ime(w), (A. 2) 

uj An 

we may rewrite the plasma frequency in terms of the dielectric function as 



OJp 



2 f°° 

— / dwwlme(w). (A. 3) 

JO 



With the help of the Kramers- Kronig relation for non-negative frequencies 

Re/M = *-V f°*J *™>, WM - r*J (A.4) 

7T J uj' 2 - uj 2 ^ Jo w - W 

we may relate Ime(u) to Ree(w) — 1. Note that for the dielectric function the polarization P(uj) is 
the response to the applied electric field E(ui), 

4ttP(uj) = (e(uj) - 1) E(uj). (A.5) 

Therefore a convenient form to write the Kramers-Kronig relation for the dielectric function is to 
include the additional —1 in the real part of e(uj), 

Ree(uj)-1 = -V duj' V 2 > = -V da/— ^ 

TT 7 uj 12 -uj 2 ituj 2 J 1-^2" 

—\ / da/a/Im e ( W ') = -^, (A.6) 

where we are interested in the high frequency regime (uj 3> V T ) since the sum rule is strictly valid 
only for uj — > oo. In experiments we deal with finite frequencies only, and thus it is possible to 
extract the plasma frequency by the following extrapolation of experimental data, 

ujI = lim (-uj 2 Re e(ui)) . (A.7) 

As explained in the main text the superconducting plasma frequency determines the frequency 
above which the superconductors becomes "transparent" in analogy with the normal metal plasma 
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frequency. The reason for this terminology follows from the fact that photons can only penetrate the 
superconductor for length scales smaller than the London penetration depth Al which corresponds 
to up s . Here, the superconducting plasma frequency should be understood with the aforemen- 
tioned analogy to normal metals in mind, described by the Drude-Sommerfeld form of the optical 
conductivity (B.6). In this case the (superconducting) plasma frequency is given by 



,2 

-2 



ne / ne 

Wp„ = 8 — / d(cj-r) -. — nr = 8 arctan(W) = Att'^^ = \7 Z . (A. 8) 

m J 1 + (wt) 2 to 



TO 



Experimentally, we cannot reach infinite frequencies and thus the sum rule is modified by a cut-off 
frequency lo c , sometimes called the partial optical sum rule of the Drude-Sommerfeld form 

ne^ Wc uj^ 

d(wr) -. — rr. — 8 arctan(o;r) = — arctan(w c r). (A.9) 

1 + (ujtY to o 4tt 

Expanding the inverse tangent function for u c — > oo, we get a series expansion in the relaxation 
rate 1 /t which reads 

arctan (w c t) = ~ — + O \ —^\ , (A.IO) 



and thus 



dwRe(7(w) rj "=f I 1- -— 1 , (A.ll) 



TT LO c T 



from which we recover (A.l) for cj c — > oo. 



B. Drude-Sommerfeld Model 

In the previous sections we made use of some facts of Drude-Sommerfeld theory, which we now 
discuss in some detail. The original Drude model is a purely classical model describing a gas of 
electrons diffusing through a fixed lattice of ions. These generate a positively charge background 
to balance the negative charge of the electron gas. Even after the discovery of quantum mechanics, 
the Drude model can still be used but some modifications are needed: The electrons are described 
by a fermionic gas obeying Fermi-Dirac statistics. This extension of the original model is usually 
called Drude-Sommerfeld model. There are several ways to arrive at the Drude formula used in 
(2.15) each with their own set of approximations. Essentially, there are three approaches: 

1. A classical approach assuming the existence of an average relaxation time describing the re- 
laxation to equilibrium after turning off an external electric field, i.e. the relaxation time 
approximation which operationally replaces the collision integral by a linearized approxima- 
tion including only the relaxation time scale and the Maxwell-Boltzmann distribution in the 
kinctic/Boltzmann transport equation. Strictly speaking, this derivation is borrowed from 
hydrodynamics where the collisions refer to a gas of weakly interacting particles. 

2. A semi-classical approach taking into account the quantum nature of the electron gas. Thus 
the free mean path is determined by the Fermi velocity of the electrons t = v-pr and only 
electrons near the Fermi surface can participate. 

3. A full microscopic approach considering an interacting electron gas employing Fermi liquid 
theory in combination with a diagrammatic expansion (and additional approximations such as 
the random phase approximation) . Here the assumptions of Fermi liquid theory, in particular 
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adiabatic continuity, are needed. Furthermore, it is useful to introduce an effective frequency 
dependent mass m* (u>) for the quasi-particles of Fermi liquid theory incorporating the effects of 
electron-phonon and electron-electron interactions as well as a frequency dependent effective 
time scale t*(u). However, in the case of impurity scattering the ratio T */m* is identical 
to the bare values since the renormalization of the lifetime cancels the renormalization of 
the mass. This can be explicitly seen in a diagrammatic calculation of the current-current 
correlator using the Kubo formula. Here the mass renormalization cancels since the diagrams 
included in the self-energy are also included in the two particle diagrams. This is in agreement 
with Fermi liquid theory stating that the current of the quasi particles is independent of the 
interaction. 

Since we do not know in general if there are really quasi-particles involved in strongly correlated 
systems (apart from heavy Fermi liquids), it seems to be a good strategy not attempt to explain 
Homes' law microscopically, but to take the "effective field theory" philosophy that focuses on the 
macroscopic rather than the microscopic degrees of freedom. In particular dealing with "Planckian 
dissipation" there is no particle transport and we are primarily interested in quantum critical 
transport at finite density. Therefore we may follow the standard derivation using the Kubo formula 
without resorting to the aforementioned technicalities. A nice exposition of this derivation can be 
found in [59]. Here we will just state the main ideas and the result. We start with an interaction 
Hamiltonian 

H int = --J-A, (B.f) 

c 

and use the Kubo formula for the conductivity, 

ct ( w < k ) = E ^ / d * < s i J (°' k ) ■ J fr k )* I s ) e ~ iwt ■ ( R2 ) 

The current density operator J can be rewritten in terms of the momentum operator. Assuming an 
exponential decay with a single time scale for all current-current correlators (relaxation time approx- 
imation) and the dipole approximation (the external electric field is constant over the characteristic 
length scale which is surely valid for vp <C c) we arrive at 

e 2 r f 

*M = — r 1 — > (B.3) 

m 1 — iuit 

where / describes the oscillator strength measuring the transition probability between two states 
given by 

f = 2 J2 \(s] Pj \s')\ 2 

Furthermore the oscillator strength can be evaluated under the assumption of free electrons i.e. 

m*ft( W ,-uv) = ^- ^d £|< s | Pj .| s ')| 2 = (p2) = ^! (B.5) 

s,s' 

since the self-interaction or electron-lattice interactions are already taken care of in the renormal- 
ization of the mass m*(u>) and the relaxation rate T~ l (ui) = ^ m p U Hty + T d-ph + T c7-cr Inserting (B.5) 
into (B.4) the oscillator strength / is given by the electron (or charge carrier) density n and thus 

np 2 T 1 to 2 1 

^M = ^ r J ^ = /^, CB.fi) 

m 1 — \UJT 47T - — 10J 

T 

using the definition of the plasma frequency (2.2) which in the limit of uj — >• reduces to (2.15). 
If we assume the correction of r e i_ p h and r e i_ i to be small compared to Ti mpur ity we can follow the 
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reasoning presented in point 3 and replace r */m" by T /m. The maximum of the real part of the 
optical conductivity (B.6) at lu = is called the Drude peak, whereas the imaginary part shows a 
maximum at w = 1 /t. 



C. Equations of Motion for the s— Wave Fluctuations 

First we expand the matter Lagrangian up to quadratic order in the fluctuations 
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as well as the Maxwell Lagrangian 
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The corresponding equations of motion for the scalar fluctuations are involved expressions so we 
will simplify them according to our needs. Since we are working in the normal phase we can set 
(f>(u) = 0. Furthermore, the only non- vanishing component of the background gauge field is A t (u) 
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with all other components being zero. In order to work out the quasi-normal modes we apply a 
Fourier transformation and assume plain wave behavior for the spatial dependence 

dwd^k 
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Thus, we end up with the following equation of motion for the scalar fluctuations 
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where h t (u), h u (u) denotes the spatial entries of the corresponding row/column and kx denotes the 
transverse vector orthogonal to the directions given by the particular equations for the component 
of a(u). For instance, we can look at the equation for one of the d — 1 spatial components, a(u) = 
a x (u)e x say, so kx = k y e y + k z e z + ■ ■ ■ and hence kx • a couples the complementary components to 
the one chosen, i.e. a y (u),a z (u), .... Furthermore, assuming k = k x e x , we see that the equations 
of the a x (u) component couples only with the at(u) component (in this case kx =0), whereas the 
equations for all other d — 2 spatial components decouple. Looking at the transverse directions e.g. 
k = k y e y + k z e z + • • • , the equation for the fluctuations in the x-direction will decouple from all 
other d — 2, but the equations for the remaining d — 2 fluctuations will couple with each other and 
a t (u). 



u 2 A' t {u) 



d-l 



u 2 uA' t (u 



3u 2 wA' t (u) . 



- 29 - 



References 

[1] P. Kovtun, D. T. Son and A. O. Starinets, Holography and Hydrodynamics: Diffusion on Stretched 
Horizons, JHEP 0310 (2003) 064 [hep-th/0309213]. 

[2] A. Buchel and J. T. Liu, Universality of the Shear Viscosity in Supergravity, Phys. Rev. Lett. 93 (2004) 
090602 [hep-th/0311175]. 

[3] P. Kovtun, D. Son and A. Starinets, Viscosity in Strongly Interacting Quantum Field Theories from 
Black Hole Physics, Phys. Rev. Lett. 94 (2005) 111601 [hep-th/0405231]. 

[4] N. Iqbal and H. Liu, Universality of the Hydrodynamic Limit in AdS/CFT and the Membrane 
Paradigm, Phys. Rev. D79 (2009) 025023 [0809.3808]. 

[5] A. Sinha and R. C. Myers, The Viscosity Bound in String Theory, Nucl.Phys. A830 (2009) 
295C-298C [0907.4798]. 

[6] J. Erdmenger, P. Kerner and H. Zeller, Non-universal shear viscosity from Einstein gravity, 
Phys. Lett. B699 (2011) 301-304 [1011.5912]. 

[7] J. Erdmenger, P. Kerner and H. Zeller, Transport in Anisotropic Superfiuids: A Holographic 
Description, JHEP 1201 (2012) 059 [1110.0007]. 

[8] A. Rebhan and D. Steineder, Violation of the Holographic Viscosity Bound in a Strongly Coupled 
Anisotropic Plasma, Phys. Rev. Lett. 108 (2012) 021601 [1110.6825]. 

[9] H. Liu, J. McGreevy and D. Vegh, Non-Fermi liquids from holography, Phys. Rev. D83 (2011) 065029 
[0903.2477]. 

[10] M. Cubrovic, J. Zaanen and K. Schalm, String Theory, Quantum Phase Transitions and the 
Emergent Fermi-Liquid, Science 325 (2009) 439-444 [0904.1993]. 

[11] N. Iqbal, H. Liu and M. Mezei, Lectures on Holographic Non-Fermi Liquids and Quantum Phase 
Transitions, 1110.3814. 

[12] S. A. Hartnoll, Lectures on Holographic Methods for Condensed Matter Physics, Class. Quant. Grav. 
26 (2009) 224002 [0903.3246]. 

[13] C. P. Herzog, Lectures on Holographic Superfluidity and Superconductivity, J. Phys. A A42 (2009) 
343001 [0904.1975]. 39 pages, 9 figures, lectures given at the Trieste Spring School on Superstring 
Theory and Related Topics. 

[14] G. T. Horowitz, Introduction to Holographic Superconductors, 1002.1722. 

[15] M. Kaminski, Flavor Superconductivity & Superfluidity, Lect. Notes Phys. 828 (2011) 349-393 
[1002.4886]. 

[16] J.-H. She, B. J. Overbosch, Y.-W. Sun, Y. Liu, K. Schalm, J. A. Mydosh and J. Zaanen, Observing 
the Origin of Superconductivity in Quantum Critical Metals, Phys. Rev. B84 (2011) 144527 
[1105.5377]. 

[17] G. T. Horowitz, J. E. Santos and D. Tong, Optical Conductivity with Holographic Lattices, JHEP 
1207 (2012) 168 [1204.0519]. 

[18] Y. Liu, K. Schalm, Y.-W. Sun and J. Zaanen, Lattice Potentials and Fermions in Holographic Non 
Fermi- Liquids: Hybridizing Local Quantum Criticality, 1205 . 5227. 

[19] S. Sachdev, Quantum Phase Transition. Cambridge University Press, 2011. 

[20] S. Sachdev and B. Keimer, Quantum Criticality, Phys. Today 64N2 (2011) 29 [1102.4628]. 

[21] J. Zaanen, Superconductivity: Why the Temperature is High, Nature 430 (Jul, 2004) 512-513. 

[22] M. Ammon, Gauge/Gravity Duality Applied to Condensed Matter Systems, Fortsch.Phys. 58 (2010) 
1123-1250. 

[23] X.-G. Wen, Quantum Field Theory of Many-Body Systems: From the Origin of Sound to an Origin 
of Light and Electrons. Oxford University Press, 2004. 



- 30 - 



[24] M. Z. Hasan and C. L. Kane, Topological Insulators, Rev.Mod.Phys 82 (Oct, 2010) 3045-3067 
[1002.3895]. 

[25] J. Zaanen, A Modern, but way too short history of the theory of superconductivity at a high 
temperature, ArXiv e-prints (Dec, 2010) [1012.5461]. 

[26] C. C. Homes, S. V. Dordevic, M. Strongin, D. A. Bonn, R. Liang, W. N. Hardy, S. Komiya, Y. Ando, 
G. Yu, N. Kaneko, X. Zhao, M. Greven, D. N. Basov and T. Timusk, A Universal Scaling Relation in 
High- Temperature Superconductors, Nature 430 (Jul, 2004) 539-541. 

[27] C. C. Homes, S. V. Dordevic, T. Valla and M. Strongin, Scaling of the Superfluid Density in 
High- Temperature Superconductors, Phys. Rev. B 72 (Oct, 2005) 134517. 

[28] A. Altland and B. Simons, Condensed Matter Field Theory. Cambridge University Press, 2006. 

[29] J. L. Tallon, J. R. Cooper, S. H. Naqib and J. W. Loram, Scaling Relation for the Superfluid Density 
of Cuprate Superconductors: Origins and Limits, Phys. Rev. B 73 (May, 2006) 180504. 

[30] D. v. d. Marel, H. J. A. Molegraaf, J. Zaanen, Z. Nussinov, F. Carbone, A. Damascelli, H. Eisaki, 
M. Greven, P. H. Kes and M. Li, Quantum Critical Behaviour in a High-Tc Superconductor, Nature 
425 (Sep, 2003) 271-274. 

[31] D. Tanner, H. Liu, M. Quijada, A. Zibold, H. Berger, R. Kelley, M. Onellion, F. Chou, D. Johnston, 
J. Rice, D. Ginsberg and J. Markert, Superfluid and Normal Fluid Density in High-Tc 
Superconductors, Physica B: Condensed Matter 244 (Jan, 1998) 1-8. 

[32] J. Mas, J. P. Shock and J. Tarrio, Sum Rules, Plasma Frequencies and Hall phenomenology in 
Holographic Plasmas, JHEP 1102 (2011) 015 [1010.5613]. 

[33] D. R. Gulotta, C. P. Herzog and M. Kaminski, Sum Rules from an Extra Dimension, JHEP 1101 
(2011) 148 [1010.4806]. 

[34] U. Giirsoy, E. Plauschinn, H. Stoof and S. Vandoren, Holography and ARPES Sum-Rules, JHEP 
1205 (2012) 018 [1112.5074]. 

[35] D. T. Son and A. O. Starinets, Hydrodynamics of R- Charged Black Holes, JHEP 0603 (2006) 052 
[hep-th/0601157]. 

[36] D. T. Son and A. O. Starinets, Viscosity, Black Holes, and Quantum Field Theory, 
Ann.Rev.Nucl.Part.Sci. 57 (2007) 95-118 [0704.0240]. 

[37] P. Kovtun and A. Ritz, Universal Conductivity and Central Charges, Phys. Rev. D78 (2008) 066009 
[0806.0110]. 

[38] S. A. Hartnoll and D. M. Hofman, Locally Critical Resistivities from Umklapp Scattering, 
Phys. Rev. Lett. 108 (2012) 241601 [1201.3917]. 

[39] N. Iizuka and K. Maeda, Towards the Lattice Effects on the Holographic Superconductor, 1207.2943. 

[40] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, Building a Holographic Superconductor, 
Phys. Rev. Lett. 101 (2008) 031601 [0803.3295]. 

[41] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, Holographic Superconductors, JHEP 0812 (2008) 
015 [0810. 1563]. 

[42] D. T. Son and A. O. Starinets, Minkowski- Space Correlators in AdS/CFT Correspondence: Recipe 
and Applications, JHEP 09 (2002) 042 [hep-th/0205051]. 

[43] M. Ammon, J. Erdmenger, V. Grass, P. Kerner and A. O'Bannon, On Holographic p-wave 
Superfluids with Back-Reaction, Phys. Lett. B686 (2010) 192-198 [0912.3515]. 

[44] T. Faulkner, G. T. Horowitz and M. M. Roberts, Holographic Quantum Criticality from Multi-Trace 
Deformations, JHEP 1104 (2011) 051 [1008.1581]. 

[45] J.-H. She and J. Zaanen, BCS superconductivity in quantum critical metals, Phys. Rev. B 80 (Nov, 
2009) 184518. 

[46] M. Ammon, J. Erdmenger, M. Kaminski and P. Kerner, Superconductivity from Gauge/Gravity 
Duality with Flavor, Phys. Lett. B680 (2009) 516-520 [0810.2316]. 



- 31 - 



[47] M. Ammon, J. Erdmenger, M. Kaminski and P. Kerner, Flavor Superconductivity from 
Gauge/Gravity Duality, JHEP 0910 (2009) 067 [0903.1864]. 

[48] J. Erdmenger, V. Grass, P. Kerner and T. H. Ngo, Holographic Superfluidity in Imbalanced Mixtures, 
JHEP 1108 (2011) 037 [1103.4145]. 

[49] F. Bigazzi, A. L. Cotrone, D. Musso, N. P. Fokeeva and D. Seminara, Unbalanced Holographic 
Superconductors and Spmtronics, JHEP 1202 (2012) 078 [1111.6601]. 

[50] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Dutta, R. Loganayagam et. al, Hydrodynamics 
from Charged Black Branes, JHEP 1101 (2011) 094 [0809.2596]. 

[51] S. S. Gubser and S. S. Pufu, The Gravity Dual of a p-Wave Superconductor, JHEP 0811 (2008) 033 
[0805.2960]. 

[52] P. Basu, J. He, A. Mukherjee and H.-H. Shieh, Hard-Gapped Holographic Superconductors, Phys.Lett. 
B689 (2010) 45-50 [0911.4999]. 

[53] S. S. Gubser, F. D. Rocha and A. Yarom, Fermion Correlators in Non-Abelian Holographic 
Superconductors, JHEP 1011 (2010) 085 [1002.4416]. 

[54] Q. Pan and B. Wang, General Holographic Superconductor Models with Backreactions, 1101.0222. 

[55] S. Sachdev, A Model of a Fermi Liquid using Gauge-Gravity Duality, Phys.Rev. D84 (2011) 066009 
[1107.5321]. 

[56] T. Nishioka, S. Ryu and T. Takayanagi, Holographic Superconductor /Insulator Transition at Zero 
Temperature, JHEP 1003 (2010) 131 [0911.0962]. 

[57] S. A. Hartnoll and E. Shaghoulian, Spectral Weight in Holographic Scaling Geometries, JHEP 1207 
(2012) 078 [1203.4236]. 

[58] M. Ammon, J. Erdmenger, M. Kaminski and A. O'Bannon, Fermionic Operator Mixing in 
Holographic p-wave Superfluids, JHEP 05 (2010) 053 [1003.1134]. 

[59] M. Dressel and G. Griiner, Electrodynamics of Solids. Cambridge University Press, 2003. 



- 32 - 



